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OPERATING  CHARACTERISTICS  OF  AN  OR-ING  AND  SELECTION 
PROCESSOR  WITH  PRE-  AND  POST-AVERAGING 

INTRODUCTION 

The  analysis  of  systems  with  or-ing  has  been  the  subject  of  considerable 
interest  [1-8].  These  analyses  have  been  in  terms  of  a  second-moment  theory 
and  have  usually  made  timely  use  of  the  Central  Limit  Theorem  at  various 
points  in  the  system,  such  as  an  accumulator  output.  A  notable  exception  to 
this  approach  is  the  exact  analysis  in  [9],  where  the  detection  and  false 
alarm  probabilities  for  a  system  with  quantizers,  or-ing,  and  an  accumulator 
were  determined  exactly.  There,  the  characteristic  function  of  the  decision 
variable  was  determined  exactly  via  a  fast  Fourier  transform  to  yield 
threshold-crossing  probabi 1 i ties. 

Here,  we  attempt  to  continue  this  theme  of  getting  exact  detection  and 
false  alarm  probabi 1  ties  for  an  or-ing  system  with  pre-  and  post-averaging. 
However,  the  analytical  situation  is  greatly  complicated  by  the  presence  of 
the  pre-  and  post-averagers:  the  pre-averager  creates  statistical  dependencies 
between  its  successive  output  time  samples  that  the  post-averager  must  operate 
on.  For  the  signal  present  hypothesis  H^,  we  thus  encounter  the  classical 
problem  of  determining  the  output  distribution  function  of  a  linear  system 
subject  to  non-Gaussian  and  non-white  noise;  since  this  is  an  unsolved 
problem,  we  must  instead  revert  to  determining  an  approximation  for  the 
characteristic  function  and/or  cumulative  distribution  function  of  the 
decision  variable.  On  the  other  hand,  for  the  signal  absent  hypothesis  HQ, 
the  characteristic  function  of  the  decision  variable  can  be  determined  in 
closed  form.  Thus  we  have  a  mixed  situation,  where  part  of  the  analysis  is 
exact  and  part  approximate. 

The  technique  utilized  here  is  not  a  second-moment  approach,  although  one 
part  of  the  analysis  does  match  moments  through  second-order,  in  order  to 
determine  the  effective  number  of  independent  samples  at  the  or-ing  input. 

The  utility  of  this  combined  approximate  analysis  is  borne  out  by  simulation 
which  corroborates  the  technique  in  a  wide  variety  of  examples  evaluated 
numerical ly. 
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The  inclusion  of  or-ing  in  a  system  in  the  usual  fashion  necessitates  an 
increase  in  the  decision  threshold,  if  the  false  alarm  probability  is  to  be 
kept  the  same  as  when  there  was  no  or-ing,  i.e.,  one  input  channel.  The 
amount  of  threshold  increase  required  is  a  function  of  the  number  of  input 
channels,  N,  and  must  therefore  be  re-calculated  every  time  N  is  changed. 

Also,  this  increased  threshold  leads  to  decreased  detection  probabilities, 
which  is  undesirable. 

Here  we  will  investigate  a  novel  or-ing  and  selection  processor  proposed* 
by  Bell  Laboratories  which  uses  a  threshold  completely  independent  of  the 
number  of  input  channels,  N,  and  yet  also  realizes  a  false  alarm  probability 
independent  of  N.  The  required  threshold  does,  of  course,  depend  on  the 
required  false  alarm  probability  and  the  statistics  of  the  noise-only  on  one 
input  channel,  in  addition  to  the  amount  of  post-averaging;  however,  it  is 
independent  of  N  and  the  amount  of  pre-averaging. 

Although  the  new  processor  utilizes  a  threshold  independent  of  N,  it 
still  suffers  a  loss  in  detection  probability  relative  to  no  or-ing,  because  a 
succession  of  decisions  must  be  reached  after  the  pre-averager  (i.e.,  in  the 
or-ing  device)  but  before  the  post-averager.  That  is,  the  nonlinear 
operations  are  not  deferred  until  after  all  the  linear  signal-enhancing 
operations  are  completed,  but  in  fact,  occur  at  an  intermediate  point  in  the 
processor.  How  well  the  new  processor  performs  in  comparison  with  a  standard 
or-ing  system  with  pre-  and  post-averaging  is  not  known;  the  latter  is 
currently  under  investigation  and  will  be  reported  on  separately. 


*U.S.  Patent  Number  4622519  has  been  assigned  to  AT&T  Bell  Laboratories, 
regarding  this  selection  processor  with  pre-  and  post-averaging. 
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DESCRIPTION  OF  OR-ING  AND  SELECTION  PROCESSOR 

The  system  of  interest  is  depicted  in  figure  1.  There  are  N  input 
channels,  and  the  input  on  the  n-th  channel  at  discrete  time  t  is  *n(t). 

The  inputs  [xn(t)}  are  stationary  and  statistically  independent  of  each 
other,  from  channel-to-channel ,  and  from  one  time  instant  to  any  other  within 
the  same  channel.  Furthermore,  the  channels  without  a  signal  in  them  are 
identically  distributed,  with  an  arbitrary  characteristic  function 
f<0)(5  ),  where  the  superscript  denotes  hypothesis  Hq,  signal  absent 
from  all  channels.  For  hypothesis  H^,  the  signal  is  present  in  one  and  only 
one  (unknown)  channel,  and  remains  in  that  channel  for  the  total  observation 
time  leading  to  a  particular  decision  at  the  system  output  w(t);  this  total 
time  is  K+M  samples  for  the  system  of  figure  1.  The  characteristic  function 
of  the  signal-bearing  input  channel  is  f^(f).  It  is  not  presumed  that 
the  signal  and  noise  are  additive  at  the  inputs  {xn(tj}  ;  all  that  is  needed 
are  the  input  characteristic  functions  f^h^  for  hypotheses  h=0  and  1  . 

(Later,  for  numerical  evaluation  purposes,  we  will  specialize  the  results  to 
Gaussian  inputs.) 

The  inputs  {xn(t)}  are  all  delayed  by  one  unit  of  time  (for  reasons  to 
be  seen  below)  and  then  subjected  to  pre-averaging  of  K  samples;  that  is,  the 
pre-averager  outputs  are 


K 

an(t)  =  xn(t-k)  for  l<n<N  .  (1) 

k=l 


Due  to  the  unit  delay,  an(t)  is  the  sum  of  the  past  K  samples  of  input 
xn(t)  on  channel  n.  The  purpose  of  this  pre-averaging  is  to  enhance  the 
signal  (if  present),  prior  to  making  a  nonlinear  decision  in  the  or-ing 
procedure. 


The  greatest-of  device  has  two  outputs.  The  first  is  the  actual  largest 
input  it  is  subjected  to  at  time  t: 


aN(tj}  =  max  {an(tjj  ~  a^(t)(t)  . 
1  <  n  <  N 


y(t)  =  max  {a^t),  ..., 


(2) 
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Figure  1.  Or-ing  and  Selection  Processor  with  Pre-  and  Post-Averaging 
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The  other  output  is  JJ(t),  as  indicated  in  (2);  it  is  the  channel  number  at 
time  t  of  the  largest  input  to  the  greatest-of  device.  Although  y(t)  is 
discarded, ^(t)  is  used  in  the  selection  procedure. 

The  channel  number  output,  Ji{ t),  is  used  to  make  a  selection  of  one  of 
the  N  inputs  at  time  t,  not  t-1.  That  is,  the  selector  output  is 


v(t)  =  x^(t)(t)  =  x(t;j£(t)), 


(3) 


which  is  one^the  current  input  samples. 

Finally,  the  post-averager  output  is  the  sum  of  M  consecutive  selector 
outputs: 


M-l 


M-l 


(4) 


This  quantity  is  called  the  decision  variable;  it  is  compared  with  a  threshold 
Afor  decisions  about  signal  absence  vs  signal  presence  in  some  input.  It  is 
important  to  observe  that  although  the  decision  variable  only  employs  M 
samples,  the  "oldest"  of  these  samples  itself  depends  on  the  previous  K  input 
samples,  by  way  of  the  pre-averaging.  Thus  if  the  processor  is  to  function 
well  under  hypothesis  H^,  the  signal  must  be  present  at  the  input  for  a 
total  of  K+M  samples,  not  just  M.  This  feature  must  be  kept  in  mind  when 
comparing  this  processor  with  others,  in  terms  of  the  detection  probability  and 
false  alarm  probabilities;  both  processors  must  be  allowed  K+M  samples  for  a 
fair  comparison. 

We  recall  from  (1)  that  the  pre-averager  outputs  at  time  t  were  unable  to 
utilize  the  current  input  samples  f x n ( t )}  ,  for  1  <_  n  <  N.  Since  the  system 
inputs  are  all  independent  at  different  times,  the  selection  procedure  will 
always  yield  a  random  variable  independent  of  those  used  in  the  or-ing.  For 
Hq,  this  means  that  post-averager  output  w(t)  in  (4)  is  the  sum  of  M 
independent,  identically  distributed,  random  variables,  each  with 
characteristic  function  f^^(f),  regardless  of  the  particular  sequence  of 
channel  values  jl(t),  j^(t-l),  . . . ,  j£(t-M+l) ,  that  occurred  at  the  selector 
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input.  Thus  the  statistics  of  w(t)  under  HQ  are  independent  of  N,  the 
number  of  input  channels,  but  depend  solely  on  f^(f)  and  M.  This 
enables  the  setting  of  a  threshold  which  realizes  a  fixed  false  alarm 
probability  ,  regardless  of  the  value  of  N. 

Under  H^,  although  the  post-averager  output  is  the  sum  of  M 
statistically  independent  random  variables?  some  of  the  selector  outputs  are 
drawn  from  the  actual  signal  channel  (if  the  pre-averager  is  functioning 
well),  but  some  are  noise-only  channel  inputs.  The  number  of  correct 
signal-channel  selections  can  be  anywhere  in  the  range  (0,M);  this  greatly 
complicates  the  analysis  for  the  detection  probability,  and  necessitates  some 
approximations. 

Behavior  of  Pre-Averager  Outputs 

Since  the  memory  of  the  pre-averagers  is  K  samples,  the  inputs  to  the 
greatest-of  device  cannot  "forget"  the  past  very  quickly.  What  this  means  is 
that  if  one  system  input  is  momentarily  very  noisy,  it  continues  to  adversely 
affect  or-ing  decisions  for  a  considerable  time.  Thus,  channel-number  output 
J?(t)  can  remain  at  the  same  value  for  several  consecutive  times,  even  if  the 
signal  is  absent.  A  sample  plot  of  this  behavior  is  depicted  in  figure  2  for 
K=10  and  N=5.  Thus,  for  example,  one  particular  channel  pre-averager  output 
dominates  from  sample  t=13  to  t=27. 

When  the  signal  is  added,  the  situation  improves,  in  that  the  signal 
channel  dominates  most  often;  see  figure  3.  However,  there  are  still  extended 
regions,  such  as  (16,26)  and  (68,90),  where  one  of  the  noise  channels 
dominates.  The  effects  of  these  extended  regions  of  dominance  must  be 
accounted  for  in  any  analysis  of  performance. 

There  are  conflicting  requirements  at  work  in  the  processor  of  figure  1. 

One  indicates  that  K  should  be  increased  so  that  good  decisions  can  be  reached 
in  the  greatest-of.  On  the  other  hand,  larger  K  can  cause  prolonged 
dominance,  as  indicated  above;  and  it  subtracts  from  the  time,  M,  available 
for  post-averaging,  if  total  observation  time  K+M  is  fixed  or  limited.  Since 
it  is  also  desirable  to  increase  M,  for  better  quality  decisions  at  the  system 

*  when  conditioned  upon  a  particular  sequence  of  ( t )  values;  see  (A-l),  (A-5), 
and  (A-6). 
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Figure  2.  Outputs  of  Pre-Averagers  for  K=10,  N=5  Channels,  Noise-Only 


Time  Sample  t 

Figure  3.  Outputs  of  Pre-Averagers  for  K=10,  N=5  Channels,  Signal  Present 
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output,  it  is  conceivable  that  there  is  an  optimum  division  between  pre-  and 
post-averaging  times,  if  K+M  is  fixed.  This  behavior  can  be  extracted  from 
the  operating  characteristics,  when  plotted  for  a  wide  variety  of  values  of  K 
and  M.  Of  course,  there  remains  the  dependency  of  any  conclusions  on  N  and 
the  input  statistics. 
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RESULTS 


Derivation  of  detection  probability  PQ  and  false  alarm  probability  Pp 
is  conducted  in  a  series  of  appendices,  under  various  degrees  of  generality. 
The  only  one  we  consider  numerically  here  is  that  of  Gaussian  inputs,  for 
which  the  input  characteristic  functions  are 


f^h)(f)  -  exp(i?y 


(5) 


-  y  3  f°r  h  =  0  and  1. 


Thus  the  two  hypotheses  have  different  means  but  the  same  variance 

2  •  * 

a  for  the  system  inputs.  We  define  the  system  input  deflection 


(6) 


The  detection  probability  and  false  alarm  probability  are  given  in  (C-19)  and 
(C-20),  where  the  various  parameters  are  also  explained.  In  the  derivations, 
it  was  necessary  to  define  the  number  of  independent  samples  at  a  pre-averager 
output  in  an  interval  of  length  M,  and  then  to  approximate  characteristic 
functions  raised  to  fractional  powers.  The  details  are  lengthy  and  are 
reserved  for  the  appendices  for  the  interested  reader. 

A  variety  of  cases  are  plotted  in  the  following  figures.  Values  of  the 
parameters  have  been  chosen  as  indicated  in  table  1  for  figures  4  -  35.  We 
have  looked  at  all  possible  combinations  of  these  choices,  for  a  total  of 
4*4*2  =  32  cases. 


N,  Number  of 
Input  Channels 


K,  Amount  of 
Pre-Averaging 


M,  Amount  of 
Post-Averaging 


2 

4 

8 

16 


1 

2 

4 

8 


32 

64 


Table  1.  Values  of  N,  K,  M  Considered  for  Figures  4-35 
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Superposed  on  most  figures  is  one  simulation  run  for  a  representati ve 
value  of  d,  equal  to  one  utilized  in  the  analysis;  for  example,  figure  4 
contains  a  simulation  result  for  d  =  .8.  These  (jagged)  simulation  curves 
closely  overlay  the  (smooth)  analysis  curves  except  in  the  region  of  rare 
events,  namely  where  the  detection  probability  is  near  unity  U-pd  «  i); 
the  random  discrepancy  in  this  latter  region  is  due  to  the  relatively  small 
and  limited  number  of  trials  employed  in  the  simulation  to  estimate  PQ, 
namely  1000  trials.  There  is  no  need  to  simulate  the  false  alarm  case  at  all, 
since  we  can  calculate  the  false  alarm  probability  exactly  in  any  region  and 
for  any  parameter  values. 

For  the  smaller  values  of  N,  the  characteristics  are  nearly  equally 
spaced;  however,  as  N  increases,  the  curves  get  more  crowded  around  the  d=0 
curve.  This  is  a  manifestation  of  small  signal  suppression,  brought  about  by 
the  nonlinear  action  of  the  or-ing  operation.  The  plots  in  figures  4-35  are 
also  generally  fairly  straight;  this  is  due  to  the  considerable  amount  of  post¬ 
averaging  used,  namely  M  =  32  or  64. 

When  M  is  decreased,  as  in  figures  36  -  39,  the  operating  characteristics 
develop  a  significant  knee,  where  the  detection  probability  temporarily 
saturates  as  the  false  alarm  probability  increases.  However,  this  undesirable 
behavior  generally  occurs  for  the  larger  values  of  false  alarm  probability, 
where  the  system  will  not  be  operated.  The  superposed  simulation  results 
confirm  the  presence  of  this  knee  in  the  operating  characteristics. 
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SUMMARY 

The  detection  and  false  alarm  probabilities  have  been  derived  for  the 
or-ing  and  selection  processor  with  pre-  and  post-averaging  and  specialized  to 
Gaussian  inputs.  Plots  of  the  operating  characteristics  for  a  variety  of 
choices  of  number  of  input  channels,  amount  of  pre-averaging,  amount  of 
post-averaging,  and  input  signal  strength  are  given.  Superposed  on  most  of 
the  analytical  results  is  a  simulation  result  for  each  case,  which 
corroborates  the  approximations  which  were  necessary  for  a  tractable 
analysis.  Other  cases  of  interest  to  the  reader  may  be  obtained  by  means  of 
the  program  listed  in  appendix  D  in  BASIC  for  the  Hewlett  Packard  9845  Desk 
Calcul ator. 
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Figure  7.  Operating  Characteristics  for  N=2,  K=2 ,  M=64 
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Figure  8.  Operating  Characteristics  for  N=2,  K=4 ,  M=32 
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Figure  9.  Operating  Characteristics  For  N=2,  K=4,  M=64 
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Figure  12.  Operating  Characteristics  for  N=4,  K=l,  M=32 
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Figure  13.  Operating  Characteristics  for  N=4,  K=l,  M=64 
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Figure  14.  Operating  Characteristics  for  N=4 ,  K=2 ,  M=32 
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Figure  18.  Operating  Characteristics  For  N=4,  K=8,  M=32 
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Figure  20.  Operating  Characteristics  for  N=8,  K=l,  M=32 
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Figure  21.  Operating  Characteristics  For  N=8,  K=l,  M=64 
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Figure  22.  Operating  Character  1 st i cs  for  N=8,  K=2,  M=32 
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Figure  23. 


Operating  Characteristics  for  N=8,  K=2,  M=64 
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Figure  24.  Operating  Characteristics  for  N=8,  K=4 ,  M=32 
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Figure  2G.  Operating  Characteristics  for  N=8 ,  K=8 ,  M=32 
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Figure  27. 
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Figure  28.  Operating  Characteristics  for  N=16,  K=l,  M=32 
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Figure  30.  Operating  Characteristics  for  N— 16,  K= 2,  M=32 
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Figure  31. 


Operating  Char acter i st i cs  for  N=16,  K=2 ,  M=64 


39 


999999 


TR  No.  6929A 


.999999 

.99999 

.9999 

.999 


.001 


.0001 


.00001 


CQ 

cn 

CD 

cn 

cn 

CD 

c n 

CD 

CD 

• 

CD 

CD 

* 

CD 

. 000001 


Figure  32. 
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Operating  Characteristics  for  N=16,  K=4 ,  M=64 
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Figure  34.  Operating  Characteristics  for  N=16,  K-8  f  M-32 
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Operating  Characteristics  for  N=4 ,  K=l,  M=1 
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APPENDIX  A.  CHARACTERISTIC  FUNCTION  OF  DECISION  VARIABLE 

With  reference  to  figure  1,  let  ^  =  {< d  ,  J.  ,  .  ..,^  )  be 

0  1  ^M-l 

the  time  sequence  of  values  of  JJ(t)  at  the  output  of  the  Greatest-Of  device 
leading  to  a  particular  value  of  the  decision  variable  w(t)  at  some  arbitrary 
time  instant.  Then  the  conditional  characteristic  function  of  w(t)  is,  since 
the  inputs  {xn]  to  the  system  are  all  statistically  independent  in  time. 


M-l 

fw(*U>  -  TTf.  (?)  , 

w  j-o 

where  the  individual  characteristic  functions  in  (A-l)  are 


(A-l ) 


(1) 


ft11!?) 

f^0)(?) 


if  signal  is  present  in  channel  no.  X- 

if  signal  is  absent  from  channel  no.  i. 

•J  J 


(A— 2) 


Recall  that  f^(f)  is  the  characteristic  function  of  a  single  channel 
input  under  hypothesis  h;  h  =  0  or  1. 


Hq :  No  signal  present  in  any  input  channel 

The  time  sequence  of  selector  outputs  {v(t)}  is  one  of  statistically 
independent  random  variables,  regardless  of  whether  sequence  {i(tj}  remains  at 
the  same  value  (channel  number)  for  several  successive  time  instants,  as  shown 
in  figures  2  and  3.  So  the  conditional  characteristic  function  of  w(t)  is, 
under  Hq5 


f(0) 


(A-3) 


independent  of  the  particular  sequenced.  Therefore  the  unconditional 
characteristic  function  of  w(t)  is 


fi0)(?)  =  [^0)(?i]M 


A-l 
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Thus,  the  cumulative  distribution  function  of  decision  variable  w(t) 
under  Hq  Can  be  easily  evaluated  [10,11]  and  the  false  alarm  probability  can 
be  accurately  calculated.  The  false  alarm  probability  is  independent  of  the 
number  of  input  channels,  N,  and  the  amount  of  pre-averaging,  K.  The  false 
alarm  probability  depends  solely  on  the  noise-only  statistics  at  a  sing.le 
channel  input,  the  amount  of  post-averaging,  and  the  threshold.  The  threshold 
need  not  be  increased  as  N  is  increased,  in  order  to  maintain  a  fixed  false 
alarm  probability.  Arbitrary  input  statistics  of  [xn(t)j  are  allowed, 
except  for  the  independent  identically-distributed  assumption. 


H^:  Signal  present  in  one  unknown  input  channel 

The  conditional  characteristic  function  of  w(t),  using  the  independence 
in  time  of  {x  (t)}  ,  is  now  given  by 


*  ( 1 ) 


<51 V  ■  [fx1*!?)]1  [fx0) 


-I  M-T 


(f> 


( A— 5 ) 


where  T  is  the  number  of  times  that  the  values  of  the  sequence^  correspond 
with  the  correct  (unknown)  signal  channel.  Thus  the  unconditional 
characteristic  function  of  the  decision  variable  w(t)  is 


M 

-  5.  Prob(T> 

T=0 

=  ]>  Prob(T)  [41}(?)]T  [fx0)(^)]|/|T  >  (A— 6) 

T=0 

where  Prob(T)  is  the  probability  of  exactly  T  correspondences  of J  with  the 
correct  signal  channel. 

Unfortunately,  evaluation  of  Prob(T)  is  very  difficult.  For  example, 
using  the  independence  of  the  channel  inputs,  the  simplest  case  of  T=M  yields 


A-2 
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Prob(T 


fo» 

=M)  = 


4-00 

r 


duM  Pa  ^ul» 


u..) 


■  u,  UM  -i 

idVl  , 

dvM  pa°)(vl’  •»  VM} 

* 

where  p^  (...)  is  the  joint  probability  density  function  of  M  adjacent 
time  samples  of  one  channel  input  to  the  Greatest-Of.  However,  probability 
density  functions  p^°^(vl5  ...,  vM)  and  p^^Up  ...,  uM)  cannot  be 
factored,  because  neighboring  time  samples  of  a  (t),  at  a  particular  pre¬ 
averager  output,  are  highly  dependent.  (If  {an(tj}  samples  were  independent 
in  time,  as  for  K=l,  this  would  simplify  to 


Prob(T=M) 


( A— 8 ) 


So  we  must  resort  to  an  approximate  analysis  under  in  order  to  determine 
Prob(T) . 


Approximate  Analysis  under 

The  plots  in  figures  2  and  3,  of  {an(t)}  vs.  time  sample  t,  reveal  that 
a  particular  pre-averager  output  can  dominate  for  many  adjacent  time  samples 
(even  if  there  is  no  signal  present).  Thus  a  typical  plot  of  J?(t)  appears  as 
in  figure  A-l. 


m 


Figure  A-l.  Sample  Plot  of  Jl{t) 

The  effective  number  of  "independent  samples",  1^  ,  of  ^(t)  in  an  interval 
of  length  M  can  be  considerably  less  than  M;  it  depends  critically  on  K,  the 
amount  of  pre-averaging.  And  since  jj(t)  depends  directly  on  the  pre-averager 
outputs  {an(t)}j|,  it  appears  that  the  effective  number  of  independent 

A-3 
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samples,  I&J  of  one  of  the  pre-averager  outputs  in  an  interval  M  should  be 
approximately  the  same  as  Ij .  This  assumption  precludes  the  dependence  of  Ig 
on  N,  the  number  of  channels;  how  good  this  presumption  is  must  be  borne  out 
by  comparison  of  analytical  results  with  simulation. 


The  effective  number,  for  a  time  sequence  an(t)  observed  for  an 
interval  M,  is  given  by  [12,  eq.  B— 12] 


I 


a 


M 


(A— 9 ) 


where  Pa(m)  is  the  normalized  covariance  of  stationary  process  an(t);  i.e., 
for  any  n. 


»a(m) 


E{an(t)  an(t-mjj 


(A-10) 


where  ap(t)  is  the  ac  component  of  an(t).  But  since  an(t)  is  the  sum  of 
K  independent  time  samples  of  xn(t),  there  follows 


1ml 


for  |m)  <  K 


»a(")  - 


0  otherwise 


( A— 1 1 ) 


regardless  of  whether  signal  is  present  or  not  in  the  channel  under 
consideration.  Substitution  of  ( A— 1 1 )  in  (A— 9)  yields 


r 


3MC 


I,  9  I  = 

Jl  a 


3KM+1-K 
3KM 

3KM+1-M2 


for  K  <  M 


for  M  <  K 


(A— 12 ) 


A  reasonable  approximation  to  (A-12)  is  1^  a  M/K  for  M  >  K;  however,  we  use 
the  more  accurate  result  (A-12)  in  the  following. 


A-4 
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Now  let  probability 

Qa  =  Prob  ^signal  channel  dominates  Greatest-Of  on  one  trial}.  (A-13) 

Using  the  statistical  independence  between  input  channels,  this  quantity  is 
given  by 


Q 


a 


P^OO 


( A— 14 ) 


this  is  the  special  case  of  (A-7)  or  (A-8)  for  T=l.  Then  in  an  interval  of 
length  M,  where  there  are  1^  effectively-independent  samples  (as  given  by 
(A-12)),  the  probability  that  the  signal  channel  dominates  the  Greatest-Of 
exactly  n  times  is 


Prob{n;  I*}  (1-Qa)  ^  "  for  0  <  n  <  1^  .  (A-15) 

We  are  presuming  1^  to  be  an  integer  for  now;  we  shall  alleviate  this 
restriction  shortly. 

Now  define 

L  =  j—  .  (A-16) 

bt 

This  is  the  average  number  of  time  samples,  in  interval  length  M,  for  which  a 
particular  channel  dominates  the  Greatest-Of.  In  terms  of  figure  A-l,  L  is 
the  average  of  lengths  L1#  |_2,  ....  We  then  approximate  the  actual  J?(t ) 
plot  in  figure  A-l  by  that  in  figure  A— 2 ;  there  are  ^  independent  steps  in 
interval  M,  each  of  length  L. 


I  ft) 


TR  No.  6929A 


Under  this  approximation,  the  random  variable  T  defined  in  (A-5)-(A-6) 
can  only  take  on  the  values 

T  =  0  or  L  or  2L  or  ...  or  M,  (A-17) 


that  is 


T  =  nL  for  0  <.  n  < 


(A— 18) 


Therefore  the  conditional  characteristic  function  (A-5)  becomes 


fi1)(fl£}  =  [fI1}^)]nL  for  0  <  n  <  h  •  (a-19) 


i  M— nL 


Employment  of  (A-15)  now  yields  the  unconditional  characteristic  function  of 
decision  variable  w(t)  under  as 


n=0 


n=0 


nL  _  M-nL 

f<0)(T)] 


=  k  [fx1}^L  +u-Qa)  [« 


(A-20) 


This  is  the  main  result  of  this  appendix. 


To  review:  the  characteristic  functions  f^^(^)  are  those  for  the 
system  inputs  {xn(tj]  under  hypotheses  h  =  0  and  1;  1^  is  given  by  (A-12) 
for  various  amounts  of  pre-averaging  K  and  post-averaging  M;  L  is  given  by 
(A— 16) ;  and  Qa  is  given  by  (A-14),  where  p^(u)  is  the  first-order 
probability  density  function  of  a  pre-averager  output  under  hypothesis  h. 
Since  an(t)  is  the  sum  of  K  independent  random  variables,  we  have 
characteristic  functions 


A-6 
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f 


(h) 

a 


for  h  =  0  and  1, 


which  can  be  used  to  determine  probability  density  functions  p^(u). 
The  fundamental  inputs  required  are  N,  K,  M,  f*  (^)»  and  f^^(f). 


( A— 2 1 ) 


For  Qa  near  1,  it  is  necessary  to  numerically  evaluate  1-Qa 
separately  from  (A-14);  via  physical  reasoning  or  an  integration  by  parts  on 
(A— 14) ,  there  follows 


1-Qa  =  (N-l)  |  du  P^1}(u)  p^0)(u)  [P(0)(u^  2  .  (A— 22 ) 

We  employ  both  (A-14)  and  (A-21)  in  the  program  that  follows. 

Special  Cases  of  Characteristic  Function  f^(f) 

By  taking  the  derivative  of  (A— 20)  with  respect  to  f  and  then  setting 
f=0,  we  find  the  mean  value  of  w(t)  under  as 

=  M  [pa  41}(i)  +  U-Qa)  x(0)(i)]  ,  (A— 23 ) 

where  x£h)(n)  is  the  n-th  cumulant  of  random  variable  r  under  hypothesis 
h.  Although  (A-23)  is  deduced  from  the  approximate  characteristic  function 
f^(£)  given  in  (A— 20) ,  this  result  is  in  fact  exact,  as  may  be  easily 
verified  by  an  independent  analysis. 

From  (A-20),  we  also  can  extract  the  variance  of  decision  variable  w(t); 
there  follows 

»‘1)(2)  ,  M  [qa  X<1J(2)  ♦  (1-Qa)  x<0)(2)  *  LQa(l-qa)  fx'^d)  -  x(°)(li[2]  . 

( A— 24 ) 

Since  the  exact  calculation  of  the  variance  of  w(t)  is  very  difficult, 
we  have  no  check  on  this  result. 


As  the  input  signal  strength  approaches  zero,  f (0)(p) ^  an(j 

X  X  5 
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P  M 

( A— 20 )  yields  (^)^.^(°)  (f)J  ,  which  checks  the  HQ  result  in 
(A— 4) . 


The  following  are  a  series  of  checks,  that  are  obviously  correct. 
If  N=l,  then  Qa=l  from  (A-14),  and* 


M 


On  the  other  hand,  if  M=l,  then  1^=1  from  (A-21),  and 


(A-25) 


:(D 


:  (  1  ) 


(0), 


'<?>  -  5a  fi  *  (1-qa>  fx 'ff)  ; 


a '  x 


(A-26) 


then  probability  density  function 


P^(u)  =  Qa  p^1}(u)  +  (1-Qa)  p^0)(u) 


giving  detection  probability 

PD  =  Prob  [w>a|Hj3  = 

-  5a  H'M  +  (Ma> 


du  p^X)(u) 


vl 


=  1-Qa  p[1)(X)  -  (1-Qa)  P^0)(A) 


(A-27) 


(A-28) 


where  A  is  the  value  of  the  threshold  in  figure  1.  Finally,  if  K=l,  then  I.  =M 
from  (A— 12 ) ,  and 


M 

=  [Qa  +  (1"Qa}  fi0)(t)]  ;  (A-29) 


this  result  and  Q  in  (A-14)  are  exact  in  this  case,  as  noted  in  (A-8). 


*  For  the  Gaussian  example  in  (C-l),  we  have  the  exact  result  for  N=l, 

r  1  x/Myxh)-A> 

Prob{w>A|Hh]  =  $(__ 
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APPENDIX  B.  NON-INTEGER  POWER  OF  A  CHARACTERISTIC  FUNCTION 

The  quantities  L  and  in  (A-20)  will  generally  not  be  integer.  So  we 
may  encounter  the  problem  that  even  if  f{f)  is  a  legal  characteristic 
function,  [f(J)]v  may  not  be,  for  rron-integer  v.  (If  it  is  a  legal 
characteristic  function,  we  ignore  the  following  procedure.)  We  circumvent 
the  illegal  characteri Stic  function  problem  by  approximating 

[f(?)]V  by  legal  characteristic  function  [f(aj)]*  f(8?)  ,  ( B— 1 ) 

where  I  is  the  integer  part  of  v,  and  a,  8  are  determined  by  equating  the 
first  two  cumulants  of  these  two  expressions  in  (B-l).  There  follows 

v  =  la  +  6,  V  =  Ia^  +  8^  .  ( B— 2 ) 

The  solution  of  these  two  equations  is 


vl 


(vl(l+l-v)) 

1(1+1) 


1/2 


3  =  v  -  Icl 


( B-3) 


Observe  that  I,  a,  8  are  functions  solely  of  v,  and  are  independent  of  the 
statistics  f(f).  We  find  generally  that  a  _>  1  and  0  <_  8  <  1;  observe  in 
figure  B-l  that  a  approaches  1  as  v  increases. 


For  v  an  integer,  there  follows  I  =  v,  a  =  1,  8=0,  as  expected.  Also, 
for  a  Gaussian  random  variable,  replacement  (B-l)  is  in  fact  exact,  as  both 
functions  yield 


12  2 

exp(ifyv  -  j  %  a  v)  for  any  v  ,  (B— 4) 

where  u  and  a ^  are  the  mean  and  variance  of  the  Gaussian  random  variable. 

We  now  apply  this  procedure  to  (A-20).  Define,  as  in  (B-3), 


B-l 


TR  No.  6929A 


us 


Figure  B-l.  Parameters  of  Approximation  (B-l) 
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1^  =  INT(L)  with  attendant  ap 
12  =  INT(I^)  with  attendant 


(B— 5) 


Then  the  interior  terms  of  (A-20)  are  replaced  by 


rf<h>(r) 


2  ^h)(Bij)  =  f{h)(?) 


(B-6) 


(B-7 ) 


The  bracketed  quantity  in  (A-20)  then  becomes 

‘  Qa  +  (1-Qa)  f(0)(D  s  f1(J)  . 

Finally,  (A-20)  becomes 


‘  [Qa  fl1)(“2t>  *  U-Qa)  f)0)(<.2j]  2  [Q3  +  <1-Qa)  f'0)(B2?n  ( B— 8) 

-  {«.  i/^w*]11  + (i-q.)  r^W1  io)^^h 

*  K[41)(«1b2E)]11  f(1)(8lB2f)  + 11-®,)  1  f(0)(BlB2?)j  (B-9) 


in  terms  of  the  input  characteristic  functions  f^.  This  form  is  useful 

f(h) 

x 

from  (B-8)  as 


when  f^‘,;  are  available  numerically.  An  alternative  form  is  available 


f'1’?)  -  [qa  fi1)(e2?)  ♦  (1-Qa)  f(0)(B2f)] 


K! 


I9-n 

Q"  (l-Qa) 


n=0 


l1)(a2^n  [fi°)(a2^12  h  •  (B-10) 


This  is  useful  if  the  products  of  characteristic  functions  can  be  handled 
easily  or  have  easily-recognized  probability  density  functions. 

B-3/B-4 
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APPENDIX  C.  GAUSSIAN  INPUT  EXAMPLE 

The  inputs  £xn(t)}  to  the  system  of  figure  1  are  Gaussian  witn 
characteristic  function 

f[h)(f)  =  exp(i^h)  -\f  a\)  for  h  =  0  and  1  .  (C-l) 

The  input  means  are  and  under  hypothesis  Hg  and 

respectively,  while  the  input  variance  is  common  value  under  both 

A 

hypotheses. 


Since  an(t)  is  the  sum  of  K  statistically  independent  random  variables, 
the  characteristic  function  of  an(t)  is 


f^h)(f)  =  exp(iIu^h)K  -  a2  K) 


(C-2) 


The  corresponding  probability  density  function  and  Jisfributien  function  are 


p(h)(u)  =  -i- 
a  ^ 


u-Ku 


(h) 


JiT, 


( C— 3 ) 


and 


(h)  T/U-Kyxh) 

p^h}  (U)  =  5*  x 


( C— 4) 


where 


u. 

0(u)  =  (2ir)-1/2  exp(-u2/2),  §(u)  =:  J  dv  0( v) 


(C-5) 


Substitution  of  (C-3)  and  (C-4)  in  (A— 14)  and  (A-22)  yields 


Qa  =  q  (VK1  d,  N-l) , 
i-Qa  =  (VTT  d,  N-l), 

where 


(C-6) 


C-I 
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and 


q(c,m)  =  j*  dx  0 ( x )  [f(x+c)]m 
q^(c,m)  s  mjdx  0(x)  [$(x)]m_1  |(x-c) 


( C— 7 ) 


(1)  (0) 


-  u. 


d  = 


(C-8) 


The  quantity  d  is  called  the  input  deflection. 
q(c,0)  =  1  and  q(c,l)  =  $(c/Y2). 


Special  cases  of  (C-7)  are 


The  characteristic  function  of  decision  variable  w(t)  under  Hq  is 
immediately  available  by  substitution  of  (C-l)  in  (A— 4 ) : 


f^0)(?)  =  exp(i^0)  M  -  \f  o\  M)  . 


Then  the  probability  density  function  of  w  under  Hg  is 

.  j_  ,r!”’ 


W  ax  WM1  ax 


and  the  false  alarm  probability  is 

00 

PF  =  Prob  Jw  >A|H01  =  J* du  p^(u)  =  $(y) 


A 


where  A  is  the  threshold  in  figure  1,  and 


m40)-a 


Y  = 


w 


( C— 9 ) 


(c-10) 


(c-ll) 


( C— 12 ) 


Linder  H,,  we  employ  (C-l)  in  (B-6)  and  get 

f‘h)<?)=  [f<h)(f)]L  =  expdf^)  U  L)  , 


( C— 13 ) 


which  is  exact  for  any  L  (as  noted  in  ( B— 4 ) ) ;  there  is  no  need  here  to  take 
the  integer  part  of  L,  as  was  done  in  (B-5)  and  (B-9)  for  the  general  case. 
C-2 


TR  No.  6929 A 


We  now  use  ( C— 13 )  in  (B-10)  and  find  characteristic  function 


n=0 


2-  /YU  „  „  *2“" 


r[Qa  exp  (1  y  m1(n)  -  \  £2  a2)  +  (1-Qa)  exp  (i*m2(n)  -  \  ?2  a2)]  ,  (C- 


where 


mi  (n) 3  M  + 


(“x1’  -  “x0))  (62  +  a2  ")  L’ 


!<">  *  “x0)  H  *  (“i1’  -  “x0>)  “2  n  L> 


2  2 
‘  *  a2M 


Here  we  have  used  (B-2)  and  ( B— 5 )  in  the  form 


14) 


( C— 15 ) 


1ft  —  i 2  ct ^  *i"  02. 

Ijj  =  1 2  “2  +  ( C — 16 ) 

and  (A-16)  in  the  form  M  =  LI^. 

The  probability  density  function  of  decision  variable  w(t)  under  is 

then 


P i1J(u) 


1  f‘2)  Oa  d-Q,/2"" 

n=0  > 


i  /u  -  mi (nA  1  /u  -  rru(n) 

*  1.7*  — r— — r- 


( C— 17 ) 


Finally  the  detection  probability  is 


C-3 
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00 


In  order  to  put  this  result  in  a  form  consistent  with  that  for  the  false  alarm 
probability  in  (C-ll),  we  use  (C-12),  (C-15),  (C-8),  and  (A-16),  to  get 


„  Io-n 

o-oa) 


* 


( C— 19 ) 


To  review:  d  is  the  input  deflection  as  given  by  (C-8);  M  is  the  amount 
of  post-averaging;  1^  is  the  effective  number  of  independent  samples  of  J!,  as 
given  by  (A-12)  in  terms  of  K  and  M;  is  the  integer  part  of  Ij  ;  a2  and 
are  the  solutions  of  ( C— 16 )  and  can  be  found  via  ( B— 3 ) ;  Qa  and  1-Qa 
are  probabilities  given  by  (C-6)  and  ( C— 7 ) ;  and  y  is  a  modified  threshold 
given  by  (C-12)  such  that  the  false  alarm  probability  is  simply 


PF  -$(Y) 


(C-20) 


As  u^)  s  then  d-»0,  and  ( C— 19 )  tends  to 

then  Q  =1  and  ( C— 19 )  reduces  to  Pp=$(Y+dyH’) . 
for  the  detection  probability  and  false  alarm 
of  this  appendix.  The  fundamental  inputs  are 


(C-20)  as  it  should.  For  N=l, 
Equations  (C-19)  and  (C-20) 
probability  are  the  end  results 
N,  K,  M,  d. 


Simulation 


Since  the  absolute  scales  of  u^,  and  ax  do  not  enter 

separately  into  the  detection  probability  and  false  alarm  probability,  but 
only  through  the  combination 


d  = 


( C— 2 1 ) 


C-4 
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we  let  B(0)  _  o  and  ax  =  1  in  a  simulation.  Then  from  ( C— 11 )  and 
(C— 12) , 


PF  =  §(-A/Yff); 


(C-22) 


this  result  is  exact,  so  there  is  no  need  to  simulate  the  false  alarm  case. 
Under  H^,  suppose  we  make  a  total  of  Ty  trials,  getting  decision  variables 


ordered  such  that  w-^  <_  w^  _< 


Wy  .  Then  the  detection  probability 


and  false  alarm  probability  are  given  by 


PD  =  Prob  [w  >A|Hj}  = 


=  1  - 


k 


PF  =  Prob  {w  >X|H0]  =  $(-w km 


at  A  =  vv^  for  1  <  k  <  Tt  .  ( C— 23 ) 


Sample  simulation  results  for  various  values  of  d 
superposed  on  the  analytical  results,  given  by  ( C— 19)  and  (C-20),  in  figures 
4  through  35. 


C-5/C-6 
Reverse  Blank 
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APPENDIX  D.  PROGRAM  FOR  OPERATING  CHARACTERISTICS 


Since  Qa  and  1-Q  are  evaluated  separately  according  to  (C-6)  and 
(C-7),  a  partial  check  on  accuracy  is  furnished  by  printing  out  their  sum. 
Also,  the  array  F(*)  holds  the  binomial  terms  in  the  upper  line  of  (C-19);  its 
sum  should  be  1  and  is  also  printed  out. 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 


T$=" Figure  39. 
N=4 
K  =  8 
M  =  8 

Dm ax =4 
D  i  n  c  = .  4 
Ds  i  rn=  1 . 2 
Tt =1000 
OUTPUT  0 


0  p  e  r  at  i  n  g  C  h  a  r  a c  t  e  r  i  st  i  c  s  '* 

NUMBER  OF  INPUT  CHANNELS 
LENGTH  OF  PRE-AVERAGER S 
LENGTH  OF  POST-AVERAGER 
MAXIMUM  VALUE  FOR  DEFLECTION 
INCREMENT  IN  d 
VALUE -OF  d  IN  SIMULATION 
TOTAL  NUMBER  OF  TRIALS 


OF 

OF 

OF 


INPUT  CHANNELS 
PRE-AVERAGERS 
POST-AVERAGER 
d 


NUMBER 
LENGTH 
LENGTH 

MAXIMUM  VALUE  FOR 
INCREMENT  IN  d 
VALUE  OF  d  IN  SIMULATION 
TOTAL  NUMBER  OF  TRIALS 


N 

K 

M 

D  m  ax 
Di  nc 
Ds  i  m 
Tt 


IN 

IK 

IM 

Drnax 
I  D  i  n  c 
;  D  s  i  m 
iTt 


OUTPUT 
OUTPUT 
OUTPUT 
OUTPUT 
OUTPUT 
OUTPUT 
OUTPUT 
INTEGER  I , J, Ns, It 
T  =  3*K*M 

IF  K< M  THEN  220 
1 1 =T/<T+1-M*M> 

GOTO  230 

I  1 =3*M*M/(T+1-K*K> 

OUTPUT  0; “EFFECTIVE  NUMBER  OF  INDEPENDENT  SAMPLES  II 
OUTPUT  0;  "  " 

1 2= I  NT ( 1 1  ) 

T= 1 2* 1 1 

A2=<T+SQR<T*< 1  +  12- 1 1 ) ) ) / < 12* < 1 2+ 1 ) ) 

B2  = I  1  - I 2*A2 


IF  N= 1  THEN  310 

OUTPUT  0;  '*  Qa","  1-Qa" Qa+O-Qa)",1' 

REDIM  F (0: I2> 

DIM  F < 0 : 64 ) 

PLOTTER  IS  " 9 8 72 A 11 
OUTPUT  705; "VS5" 

LIMIT  10,180,0,235  !  LABELLING 

MOVE  51,25 
DIM  T$ [ 37 ] 

N$  = "  for  N= 11 


K$= " ,  K= " 
M$= " ,  M=M 
B$=CHR$<8> 
CSIZE  3, . 5 
LORG  5 


LABEL  T$;N*; B$;n;b*;k$;B*; K; b$;m$; B$;M 

LIMIT  10,180,65,235  !  GRID 

LOCATE  12,98, 12,98 

DATA  1 E-6 , IE-5,  IE-4, .001 , . 01 , .  1 , . 5 

DATA  . 9, . 99, . 999, . 9999, . 99999, . 999999 

DIM  P  <1 :  1 3 ) , Y  d : 1 3  > 

READ  P<*> 


_  H 


I  1 


SUM  F<*> 


D-\ 
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518 

520 

530 

548 

550 

560 

570 

5S0 

590 

600 

610 

620 

630 

640 

650 

660 

670 

630 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

310 

820 

330 

840 

850 

360 

370 

380 

390 

900 

910 

920 

938 

940 

950 

960 

970 

980 

990 

1000 


!  NORMAL  ?RO£A£>tLirr 


SCALES 


FOR  1=1  TO  13 
YU  >  =  F  N I  n y p h i  CPC  I  >  ) 
NEXT  I 
S=Y  < 1 > 

B  =  Y  < 13) 

SCRLE  S ,  B ,  S ,  B 
FOR  1=1  TO  13 
MOVE  S  ,  Y  C  I ) 

DRAW  B,YCI> 

NEXT  I 

FOR  1=1  TO  13 
MOVE  YCI),S 
DRAW  Y< I ) , B 
NEXT  I 
MOVE  8,3 
DRRW  B , E 
CSIZE  3, . 5 
LORG  1 
MOVE  . 4, . 2 
LABEL  Hd=0" 

LORG  8 

FOR  1=1  TO  13 
MOVE  S+. 1, YCI) 

LABEL  PCI) 

NEXT  I 

R  = . 5*<YC7)+Y<8) > 
MOVE  S- 1 ,  R 
CSIZE  4, . 5 
LORG  3 
LABEL  "P" 

MOVE  S-l , R- . 1 
CSIZE  2, . 5 
LORG  6 
LABEL  MD" 

CSIZE  3, . 5 
LDIR  PI/2 
LORG  8 

FOR  1=1  TO  13 
MOVE  Y  < I > , S  + . 0 5 
LABEL  PCI) 

NEXT  I 
MOVE  R ,  S-  1 
CSIZE  4, . 5 
LDIR  0 
LORG  8 
LABEL  "P" 

MOVE  R, S-l . 1 
CSIZE  2, .5 
LORG  6 
LABEL  "F" 
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1010  Sk  =SQR  <  K  >  !  ANALYTICAL  RESULTS 

1020  Sm=SQR<M> 

1030  P 1 3=P (13) 

1040  RED  I M  Xp<l:Tt),Yp<l:Tt> 

1050  DIM  Xp<l: 1000), Yp<l: 1000) 

1060  FOR  D s  =  D l n c  TO  Dm ax  STEP  Dine  !  DEFLECTION  d  IN  (6) 

1070  Dk=Ds*Sk 
1080  D  m  =  D  s  *  S  rn 
1090  IF  N > 1  THEN  1130 
1100  MOVE  S , S  +  Dm 
1110  DR AW  B,E  +  Srn 
1120  GOTO  1450 

1130  Qa=FNQs  < Dk , N- 1 >  !  <C-6> 

1140  Qal=FNQsl <Dk , N-l)  !  CC-6) 

1150  F  <  0  >  =Qal A 1 2 
1160  FOR  Ns=  1  TO  12 

1170  F(Ns)=F<Ns-l>*< I2+l-Ns)*Qa/<Ns*Qal )  !  EINOMIRL  TERMS  IN  CC-19) 

1180  NEXT  Ns 

1190  OUTPUT  0; Qa, Qal , Qa+Qal , SUM<F> 

1200  T=Dm/Il 
1210  Ta2=T*A2 
1220  Tb2=T*B2 
1230  J=0 

1240  FOR  Gam  =  S  TO  E  STEP  <E-S>-'50 

1250  J=J+1 

1260  X  p  <  J )  =  G  am 

1270  G  2  =  G  am  +  T  b  2 

1280  Pd=0 

1290  FOR  Ns=0  TO  12 
1300  V  =  F  N  P  h i <G2+Ta2*Ns) 

1310  V  1  =  F  N  P  h  i <Gam+Ta2*Ns> 

1320  P d  =  P d  +  F ( N s ) * <  Q a* V  +  Q a 1  * V 1 >  I  <C-19) 

1330  NEXT  Ns 
1340  D  =  Pd-  -  5 

1350  IF  flBSCDX . 49999999999  THEN  1380 
1360  Yp< J)=100*SGN(D> 

1370  GOTO  1390 

1380  Y  p  <  J  )  =  F  N  I  n  u  p  h  i (Pd) 

1390  IF  P d  > P 1 3  THEN  1410 
1400  NEXT  Gam 
1410  FOR  1=1  TO  J 
1420  PLOT  Xp < I > , Yp < I > 

1430  NEXT  I 
1440  PENUP 
1450  NEXT  Ds 
1460  PENUP 


D-3 


^  No. 

1470 

1430 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1  650 

1660 

1670 

1630 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1730 

1790 

1300 

1310 

1820 

1330 

1840 

1850 

I860 

1870 

1380 

1390 

1900 

1910 

1920 

1930 

1940 

1950 

I960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 


6^2?  A 


L4=-L0G  <  4 ) 

R=P 1/130 
Km=K+M 

RED IM  X<  1 : Km, 1 : N> , M< 1 : Tt > 

Din  X < 1 : 72 , 1 : 1 6 > , N < 1 : 1 000  > 

FOR  It  =  l  TO  Tt 

FOR  1=1  TO  Km 

FOR  Ns= 1  TO  N  STEP  2 

H=1E-5*INT<1E5*R> 

L  =  R-H 

U  =  FRfiCT  <31622 U0 000 0*L> 

V  =  F  R  A  C  T  <777021*H> 
W=FRACT(777021*L) 

R  =  FRFICT  (U+V+W+1E-1 1  ) 

V  1  =  R  -  .  5 

H=1E-5*INT(1E5*R) 

L  =  R-H 

IJ  =  FRRCT  <31622000 0O0*L> 
V=FRACT<7 77021 *H> 
W=FRACT(777021*L> 

R=FRRCT (U+V+W+1E-1 1  > 

V  2  =  R  -  .  5 

S2=V1*V1+V2*V2 
IF  S2  > . 25  THEN  1550 
Q2=(L4-L0G(S2>)/S2 
Q2=SQR ( Q2+Q2 ) 

X ( I , Ns  >  =V 1 *Q2 
IF  Ns+ 1 >N  THEN  1760 
X< I , N  s  + 1 )  =  V  2  *  Q  2 
NEXT  Ns 

X< I , 1 >=X< I , 1 )+Dsi m 

NEXT  I 

Wt=0 

FOR  I =K+ 1  TO  Km 
Vt  =-9E99 
FOR  Ns= 1  TO  N 
R  =  0 

FOR  J=1  TO  K 
R=R+X< I-J, Ns) 

NEXT  J 

IF  Yt  >fi  THEN  1900 

Yt  =fl  ! 

Lt =Ns  !  CHANNEL  NUMBER  OF 

NEXT  Ns 

Wt=Wt+X< I , Lt > 

NEXT  I 
W( It )=Wt 
NEXT  It 
HAT  SORT  W 
F=l/Sm 

FOR  1 1  =  1  TO  Tt 
Pd=l-( It-. 5)/Tt 
Yp(It)=FNInvphi (Pd) 

Xp( It )=T=-W( It  >*F 
IF  T<  S  THEN  2030 
NEXT  It 

FOR  1=1  TO  HIN(It,Tt> 

PLOT  X p ( I ) , Y p ( I ) 

NEXT  I 
PENUP 
END 
1 


!  SIMULATION  RESULTS 
!  SEED  FOR  RANDOM  NUMBERS 


I  GENERATE  TWO 
!  INDEPENDENT 
!  NORMALIZED 
!  GAUSSIAN 
1  RANDOM  VARIABLES 
!  VIA 

!  ACCEPTANCE- 
!  REJECTION 
!  METHOD 

!  SIGNAL  IN  CHANNEL  1 


!  PRE-AVERAGER  OUTPUT 

!  OR-ING  OPERATION 
LARGEST  PRE-AVERAGER  OUTPUT 
LARGEST  PRE-AVERAGER  OUTPUT 

!  POST-AVERAGER  OUTPUT 


* 
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r 


2090  DEF  FNQs  <  Cs , Ms  >  !  <C-7>,  LINE  1 

2100  INTEGER  K 

2110  IF  Ms  =  0  THEN  RETURN  1 

2128  IF  Ms  > 1  THEN  2150 

2130  V  =  F  N  P  h i < Cs* . 707 1 0678 1 1 87 > 

2140  RETURN  V 
2150  fla-7 
2160  B=7 

2170  S=  <  FNS <  fl , Cs , Ms ) +FNS  <  B , Cs , Ms  >  >  * . 5 

2130  N  =  2 

2190  H=  <  B- A ) * . 5 

2200  F=<B-fi>/3*. 398942280401 

2210  T  =  0 

2220  FOR  K= 1  TO  N-l  STEP  2 
2230  T=T+FNS<fl+H*K, Cs, Ms) 

2240  NEXT  K 
2250  S  =  S  +  T 
2260  Vo= V 
2270  V=  <  S  +  T  >  *F 
2280  PRINT  V 

2290  IF  RBSa-Vo/VX  =  lE-9  THEN  RETURN  V 

2300  N  =  N  +  N 

2310  H=H* . 5 

2320  F  =  F* .  5 

2330  GOTO  2210 

2340  FNEND 

2350  ! 

2360  DEF  FNQs 1 <  Cs , Ms)  !  <C-7>,  LINE  2 

2370  INTEGER  K 

2380  IF  Ms=0  THEN  RETURN  0 

2390  IF  Ms  > 1  THEN  2420 

2400  V  =  F  N  P  h i < -Cs* . 707 1 06781 1 87 > 

2410  RETURN  V 
2420  fi=-7 

2430  B=7 

2440  S=<FNS1  <fi,  Cs,  Ms>+FNS1  <B,  Cs,  Ms)  >*.5 

2450  N=2 

2460  H=  <  B-fl )  * .  5 

2470  F= < B-fl X3*Ms*. 398942280401 
2480  T  =  0 

2490  FOR  K= 1  TO  N-l  STEP  2 
2500  T=T+FNSl<fi+H*K, Cs, Ms) 

2510  NEXT  K 
2520  S=S+T 
2530  Vo=V 
2540  V= ( S+T ) *F 
2550  PRINT  V 
2560  Dv=flBS<V-Vo> 

2570  IF  <  Dv<  =  lE- 10)  AND  C Dy/V<  =  1 E-6 >  THEN  RETURN  V 

2580  N  =  N+N 

2590  H  =  H* . 5 

2600  F=F* . 5 

2610  GOTO  2480 

2620  FNEND 

2630  ! 


D-5 
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2640 

2650 

2660 

2678 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2320 

2330 

2840 

2350 

2360 

2370 

2830 

2390 

2900 

2910 

2920 

2930 

2948 

2950 

2960 

2970 

2930 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 


DEF  F N S <  X , C s , M s >  t  INTEGRA N D  OF  <C- 

T  =  F  N  P  h  i  <  X  +  C  s  > 

RETURN  E  X P ( - . 5  *  X  * X ) * T  M  s 
FNENB 

j 

DEF  FNS1 (X, Cs, Ms)  !  INTEGRAND  OF  <C- 

T=FNPh i <X> 

T 1 =FNPh l <  X-Cs ) 

RETURN  E X P < - . 5*X*X>  *TM  Ms- 1 )*T1 
FNEND 

BEFFNPhi<X)  i  <  C - 5  > 

INTEGER  J 

IF  ABS <  X  >  >5 . 14  THEN  3030 
A=. 28284271 2475*X 
C=COS<A> 

S  =  S I  N  <  A  ) 

B  =  C  +  C 
A=B*C- 1 

C=A*< 1 . 2536751 E-18+B+7. 1 0005E-20+A*7 . 451 7E-21 ) 
C=A*  < 1 . 533 4234 2 5 E-l 6  +  B* 1 . 01 64 9277 E- 1 7+0 
C  =  A  *  < 1 . 36760444757E-14+B+1 . 0601364636E-15+O 
C  =  A  *  <  8  . 39786526722E-13+B+8. 06068333945E-14+C} 
C=A*<4. 226 161443 18E-1 1+B*4. 46968229249E- 1  2  +  0 
C  =  A*< 1 . 466606 1 4234 E-9+B* 1 . 308485378 1 @E- 1 0+C > 

C = A  *  <  3 . 7225234 9369E-3+B*5 . 34275027603E-9+C > 

C = A  *  <  6  ■ 91927520325E-7+B*l . 1 5330990944E-7+C > 

C  =  A  *  <  9 . 43231 169833E-6+B*l . 320663 1 6364E-6+0 
C=A*<9. 44909263810E-5+B*2. 1 0404583073E-5+C > 

C  =  A*< 6. 971 8379240 3E-4+B+1 . 732280 16255E-4+C) 

C=A*<3. 301 5076793 5E-3+B*l . 1 0860645342E-3+O 
C=R*< . 0 1 53985726 157+B*. 0050790696 1220+C) 

C=A*< . 0 467 75523 4325+ B*. 0 1 72439625837+C ) 

C=A*< . 108630 245023+ B*. 043977338 1 94 1 +C ) 

C  =  A  *  < . 20 133 97472 65+B*. 0869894549959+O 
C=A*(.330501521917+B*. 144227226362+0 
C=. 7 0322 500 2 7 44 +B*. 24 7255 1631 40+C 
Phi =. 5+. 04501 58 158079*X+. 5*S*C 
RETURN  Phi 

IF  X  >7  THEN  RETURN  1 

N  =  M A X <  6 ,  INT(69/ABS<X>  > ,  I  NT ( 525/ < X*X ) ) )  +  l 
S  =  A=  1 
B=l/X 
C  =  B*B 

FOR  J=1  TO  N 
A=< 1-J-J)*A*C 
S  =  S  + A 
NEXT  J 

Phi  =.  3  9  3  9  4  2  2  3  0  4  0 1  *  E  X  P  < - . 5  *  X  *  X  >  *  A  B  3  <  B  >  *  S 
IF  X>0  THEN  Phi = 1 - P h i 
RETURN  Phi 
FNEND 


7),  LINE  1 


7 > ,  LINE  2 


D-6 


3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3330 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3430 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 
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DEF  FNInvphiCZ)  !  INVERSE  FUNCTION  TO  PHI  IN  <C-5) 

INTEGER  N 

DIM  T<0: 28) , R<@: 20 > 

DATA  . 992885376619, . 120467516143, . 0 1 6078 1 99342 1 , . 002686704437 1 6 
DATA  .  49963473024E-3,  .  988982 1 86E-4 ,  .  2039 1 8 1 276E-4 ,  .  432727 162E-5 
DATA  . 9380814  IE-6, . 20673472E-6, . 461597E-7, . 1041663E-7, . 23715E-8 
DATA  . 54393E-9, . 12555E-9, . 2914E-10, . 679E-1 1 , . 159E-1 1 , . 37E-12 
DATA  . 912158803418, 0 1 626628 1 8677 , . 43355647295E-3 , . 2 1 443857007E-3 
DATA  . 262575 108E-5, - . 302 1 09 1 05E-5 , - . 1240606E-7, . 6240661E-7,  54912E- 

DATA  1 4232 1 E-8 , . 3438E-10, . 3358E-10, 146E-11 , 31E-12, . 5E-13, . 2E-1 

DATA  . 956679709020, 0231070043091 , 00437423609751 , -■ 576503 42 26 5 E -3 
DATA  1 0961 0223 IE-4, . 25 1 0854702E-4 , . 1 056233607E-4 , . 27544 1233E-5 
DATA  . 43248450E-6, 2053034E-7, 43891 54E-7, 1 76340 1 E-7 399 1 29E-8 
DATA  1 3693E-9 ,  .  27292E-9,  .'13282E-9,  .  3133E-10,  .  167E-1 1 ,  204E-1 1 

DATA  965E- 1 2 , 22E-12 
X=Z+Z-1 
A=B=ABS<X> 

IF  A >= . 8  THEN  B  =  SQR C -LOG < 4*Z* < 1 -Z ) ) > 

IF  AC. 8  THEN  3450 
IF  AC. 9975  THEN  3410 
Nniax  =  20 
RESTORE  3270 

Y=-. 55945763 1330*E+2. 28791571626 
GOTO  3480 
Nrnax  =  1 5 
RESTORE  3240 

Y=-l . 54381 384237*B+2 . 56549012315 

GOTO  3430 

N  rn  ax  =  1 8 

RESTORE  3200 

Y=12. 5*Z*<2-1 >+2. 125 

REDIM  A ( 0 : N m ax ) 

READ  A  <  * ) 

Y2  =  Y  +  Y- 
T  <  0  )  =  1 
T  <  1  >  =  Y 

FOR  N=2  TO  Nmax 
T (N)=Y2*T (N-l >-T  <  N  -  2  ) 

NEXT  N 
R  =  0 

FOR  N-Nmax  TO  0  STEP  -1 
R  =  R  +  A  ( N  )  *T  <  N  > 

NEXT  N 

In vphi=SGN<X>*B*R*l. 41421356237 

RETURN  Invphi 

FNEND 


d-t/d-  y 

'Reverse 
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